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Abstract 

We calculate the exact zeros of the partition function for a continuum system 
where the probability distribution for the order parameter is given by two 
asymmetric Gaussian peaks. When the positions of the two peaks coincide, 
the two separate loci of zeros which used to give first-order transition touch 
each other, with density of zeros vanishing at the contact point on the positive 
real axis. Instead of the second-order transition of Ehrenfast classification as 
one might naively expect, one finds a critical behavior in this limit. 

PACS numbers: 64.60.Fr, 02.30Dk, 02.50Cw, 05.70.Fh 



Typeset using REVTgK 



* j ul@hep 1 . c . u-t okyo .ac.jp 
^ kclee @ phya . snu . ac . kr 



I. INTRODUCTION 



It has been a central theme since the discovery of statistical mechanics to understand 
how the analytic partition function for a finite-size system acquires a singularity in the 
thermodynamic limit if the system undergoes a phase transition [0]. The Lee- Yang theory 
0] has partly furnished the answer to this quest. They proposed a scenario where the zeros 
of the partition function form a line and cut across the real axis. They showed that the 
discontinuity in the 1st order derivative of partition function is proportional to the angular 
density of zeros, using an analogy with the two dimensional electrostatics. Then they proved 
this scenario for Ising-like discrete systems under very general conditions. They could show 
that the zeros were distributed on a unit circle in this case. 

There have been many attempts to generalize the "Lee- Yang circle theorem" ever since. 
Fisher initiated a study of zeros of the partition function in the complex temperature 
plane and extensive studies of zeros of the partition function in complex temperature plane 
followed In these works authors considered continuous phase transitions or critical 

points. 



The conceptual basis of the Lee- Yang circle theorem was finally clarified in ref. |T3 
by considering the first-order transition of a system with more general continuous degrees 
of freedom, with a doubly peaked probability distribution for the order parameter. Since 
the Ising-like models considered by Lee and Yang would be described by two symmetric 
Gaussian peaks in the thermodynamic limit, this result provides a simple conceptual basis 
for Lee- Yang unit circle theorem. Furthermore it is a generalization since general asymmetric 
configurations were considered, whose zeros form a curve which is not a unit circle in general. 

One interesting problem to consider is what happens when the positions of the two 
Gaussian peaks coincide. Since this is the limit where the latent heat / vanishes, one might 
naively expect that the system would exhibit a second-order transition of the Ehrenfest 



classification |TJ], where there is a finite discontinuity in specific heat but no latent heat. 
(Fig-S). 

However, when we consider the exact zeros of the partition function for the system with 
two Gaussian peaks, we find there is a branch of zeros other than the one described in 



ref. |T^. For / 7^ 0, this branch can be neglected, since for generic systems the Gaussian 



approximation breaks down at this point due to the contributions from the higher order 
cumulants. However, for / = this is no more true and we have to take this branch into 
account. Because of this, the system exhibits a critical behavior instead of the second-order 
transition. 
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FIG. 1. The energy density as a function of the reduced temperature at the second-order 
transition. The branch (a) is for Ac > 0, and (b) is for Ac < 0. 

II. LOCUS AND DENSITY OF ZEROS 

We consider a canonical partition function of a continuum system, which can be written 
as in ref . [[T^ 

/oo 
e'^f{x)dx (2.1) 
-oo 

where the probability density function is given by 

/oo 
Q{x/(3o)e-''dx, (2.2) 
-oo 

t = 1 — P/Pq, X = PqE, Q{E) is the density of states at energy E and Pq is the inverse 
of the transition temperature we are interested in. When one is interested in field driven 
phase transition, one may replace the energy E by the magnetization M and the inverse 
temperature /5 by the magnetic field H in case of magnetic systems, so on. We investigate 
the locus of zeros on the plane of complex temperature z = re*^ = e*. 

Now consider the case of interest in this paper, when f{x) is given by summation of two 
Gaussian peaks up to normalization, 

/(^) = -T^e^Ph . 2 ]+^^^^M- S \ (2-3) 
v27r(Ti 2crf v2vr(T2 ^^2 

These two peaks represent two different phases of the system. When /xi ^ \i2 the system 
undergoes the first-order transition. Let us denote two Gaussian functions by f\{x) and 
/2(x). Since we can relabel f\ and /2, and redefine a, we may assume m = (/i2 — > 

without loss of generality. We then have 

/•oo 

M{t)= / e*"[/i(x) + a/2(x)]rfx 

J oo 



exp(V'i) + exp(?/'2 



(2.4) 



where 



exp{tpi{t)) = J e*"^ fi{x)dx 
exp(V^2(t)) = / e^''af2{x)dx. 



(2.5) 



From the expressions above, one can easily see that the locations of zeros are given by the 
solutions to the following equation as in ref. 



Mtk) - Mtk) = 2ilk = i{2k + 1)71 



(2.6) 



where k runs through all the integers. For the double Gaussian distribution the equation 
above can be rewritten as 



1 



ilk = - \na + m{tk) + —{t 



2 ^ 2 

where = (cr| — crf)/2. This equation is quadratic and easily solved. The solutions are 



(2.7) 



4 



m . \/2 / , ^ 1 , m? 



^ ± -; — r x Hlk In a + 

a V 2 2^2 



^2 



m \h\ ... ,j ^|Afc, 
-:r^ ± -p— ± ^ sign(4)| — 



(2.8) 



where 



In a. 



' + 11 



In a. 



(2.9) 



Note that there are two branches of solutions. One passes through the transition point t = 
in the thermodynamic limit and the other does not, so the latter was implicitly discarded in 
ref. ITBIl. As we will see, the second branch closes in toward t = as we take the limit / 0. 



Now we redefine the variables 



m 



a' 



m 

Nyk 
2 



(2.10) 



and consider the thermodynamic limit N oo. We then get 



ln(r,) = ^{h 



m ^ 
a" 



Ok = ^{tk) = ±sign(4 



a 



±sign(?/^ 



1. ^ ^4+ yl 



\ ^ 4^Ac^ (Ac)2 2(Ac)2 



(2.11) 



The terms involving In a are finite size corrections and vanish in this hmit. We solve the 
second equation of ( |2.11| ) in terms of yk to get 



yk = ±ekhl + 9l 



,Ac 



(2.12) 



We substitute (|2.12| ) into the first equation of (|2.11| ) to get the locus of zeroes, 



r± = exp 



± 



Ac lAcI 



(2.13) 



We can also obtain the angular density of zeros. By taking formal derivative with respect 
to the integer k, we get 

d9k 



dk 



^1/ \l 4{ac) ^ (Ac)2 2(Ac)2 



X ■ 



_ 27r/^l + e2(Ac)2 

~ N[P + 2(Ac)2^2] 
Therefore, the angular densities of zeros g± of two branches are given by 



2tt 

2ng^{e) ^ - 



dk 



de 



^1 + 2(^)2^2 

1 + (4^)2^2 



(2.14) 



(2.15) 



III. FIRST-ORDER TRANSITIONS 

We will now consider both loci of zeros of the partition function at first-order transition. 
Note that all the quantities above depend only on the ratio //Ac except for a overall factor 
of / in front of 5'(6').Q When //Ac ^ 0, we get the first-order transition. This is the case 



^When both / and Ac are zero these quantities are ill-defined and we can no longer use the 
Gaussian approximation. One then has to take into account higher order cumulants. 



considered in ref. [|1^. There only the locus of zeros near the transition point t = were 
treated carefully since these were the only things of interest. In fact, for generic systems we 
expect that the Gaussian approximation breaks down away from the transition point t = 
due to the higher-order cumulants. 

Let us elaborate on this point. The locus of zeros cross the real axis at t = and 
t = —21/ Ac, indicating there are two phase transitions. This can be easily understood. The 
probability density at arbitrary temperature is given by 

tx -{x-xif/{2af) tx -{x-X2Yl{2ul) 
= — U — 

-{x-xi-aftf /{2al)+xit+alt'^ /2 ^-{x-X2-a^t)^ / {2(7^)+X2t+ajt'^ /2 



We see that for nonzero t the positions of the peaks are shifted, and also the relative weights 
change. We see that the position of the peak for large ai gets shifted by a larger amount 
for given temperature change, consistent with the fact that it has larger specific heat. The 
weight of the peak 1 relative to the peak 2 is given by: 

C^2 _ 2\ 

wl/w2 = exp{^ + {xi - X2)t) (3.2) 

By construction, at t = 0, the weight of two Gaussian peak is equal. Assuming a2 > cti, we 
see that for t > the peak labeled by 2 dominates. When t becomes slightly negative, then 
the peak 1 dominates. Also the positions of the Gaussian peaks get shifted to left, but the 
peak 2 moves faster. For t < —{xi — X2)/{<jI — 0^2) the peak 2 goes to the left of the peak 
1. At t = —21/ Ac, the weight of the peak 2 become equal to that of 1 again, and the peak 
2 is dominant for t < —21/ Ac. Therefore at this temperature there is another first-order 
transition with latent heat I and specific heat change —Ac. We can make similar arguments 
for (Ti > o"2. This process is depicted in Fig.^ 

This mechanism works only if we trust that the Gaussian form given in ( p.3| ) is exact. 
However, for a generic system, this is just a leading truncation of the cumulant expansion 

exp[-Nf{x)] = exp[-Ar(/(xo) + ^^{Ax)' + ^^{Axf + ■■■)], (3.3) 



so the higher order cumulants can be ignored only when Ax << 0(1/ viV). But at the first- 
order transition at t = —21/ Ac, the system is dominated by the peaks which are located 
at the distances of 0(1) from the positions of the peaks at t = 0. Therefore the higher- 
order cumulants would contribute, and we cannot trust the picture above. However, when 
Ax << 0{1/ \/N), or when the higher order cumulants are extremely small due to some 
reason, transition at t = —21/ Ac cannot be neglected anymore. In particular, in the limit 
/ — >■ 0, the second branch touches the first branch, preventing the system from exhibiting 
the second-order transition. 
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FIG. 2. Qualitive behavior of e*^f](x) versus x, where we take ^1^2 = i5 in this example, (a) 
For t > the peak 2 dominates (b) At t = the weight of the two Gaussian peak is equal by 
construction, meaning the areas under the curves are same, (c) As t decreases below 0, the weight 
of peak 1 becomes larger. The first-order transition has occurred with the latent heat I and the 
discontinuity in the specific heat Ac. The positions of the peaks begin to get shifted to th left, with 
peak 2 moving faster, (d) For t < —21/ Ac, the weight of the peak 2 becomes larger than that of 1 
again. At this point the peak 2 is at the left of the peak 1, so this is another first-order transition, 
with the latent heat I and the discontinuity in the specific heat —Ac. (e) Schematic diagram of 
energy versus reduced temperature for the system with double Gaussian peaks. Note that there 
are two first-order transition with the same latent heat but opposite sign for the discontinuity in 
the specific heat. The transition at t = —21 /Ac is discarded for generic systems. 



The behaviors of the loci of zeros for various values of Ac// are depicted in Fig. 1,2 and 
3 in complex z = exp(t) plane. t+ is the outer curve and t- is the inner curve. Only the 
zeros in the first Riemann sheet are shown. When Ac/l > (< 0), t+ passes through 
t — 0, and becomes unit circle as one approaches the symmetric limit, Ac/l — > 0. This is 
consistent with Lee- Yang's unit circle theorem. The other branch t_ degenerates to 
origin. (goes to infinity.) The loci intersect the real axis orthogonally as long as / 7^ 0. 
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FIG. 3. (a) Ac/Z = 1. As Ac/Z 0, the outer curve becomes a unit circle and the inner curve 

degenerates to origin. 

(b) Ac/l = —1. As Ac/l —> 0, the inner curve becomes a unit circle and the outer curve goes away 
to infinity. 
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FIG. 4. (a) cr^/m = Ac/Z = 10. The zeros for finite N are also plotted as dots, for 
2(7^ = N/S.C = 10 and a = 1.0. The dots deviate from the curve if a 7^ 1.0. 

(b) Magnification of the box in (a). The horizontal line indicates the real axis, and the vertical 
line is given by Re{z) = 1. Note that the outer curve passes through the point z = 1.0, and the 
angles between both loci and the real axis are 90 degrees. 
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FIG. 5. (a) Ac// = -10. The dots indicate zeros for 2a^ = NAc = -10 and a = 1.0. 
(b) Magnification of the box in (a). It is the inner curve which passes through z = 1.0 in this case. 

IV. L ^ LIMIT AND THE CRITICAL BEHAVIOR 

The limit //Ac = may be considered as the opposite limit from the symmetric case 
Ac// = 0. Now the two loci t± which were separate when / ^ 0, touch each other at ^ = 
and form a single curve. (Fig.P). Their loci are given by 

r± = exp{±\0\) (4.1) 

The density of zeros is 



2ng{9) = 2n{g+{9)+g_{9)) = 4Ac\9\. 



(4.2) 
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FIG. 6. (a) / = 0. The zeros for A^|Ac| = 10 are plotted, with the same notation as before, 
(b) Magnification of the box in (a). Note that the two loci are now joined to form a single curve, 
and intersects the real axis at z = 1.0 with the angle of 45 degrees. Both Ac > (Fig.^ and 
Ac < (Fig.|5|) approach this limit as / — > 0. 

Note that g{6) is zero at = 0, consistent with the fact that the first derivative of the 
partition function has no discontinuity. The loci intersect the real axis with the angles of 
45 degrees, and the intermediate region dominated by peak 1 with the smaller specific heat, 
which used to separate two domains dominated by peak 2 with the larger specific heat, 
touches the real axis at just one point. Therefore the system is dominated by the peak 2 
for t 7^ 0, and the peak 1 has same weight as the peak 2 only at t = 0, when their positions 
coincide. The qualitative behaviors of the two peaks for t > and t < are same as the 
ones depicted in Fig. 2(a) and Fig. 2(d). 

Therefore the system is exhibiting a critical behavior where it is just on the verge of 
making a phase transition. However, in contrast to many familiar examples of critical 
behavior, the specific heat near t = remains finite instead of blowing up. At this stage it 
is not yet clear whether there is an example of a discrete system whose critical behavior at 
the thermodynamic limit can be described by this model. 
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